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Within mean-field theory we calculate the scaling functions associated with critical Casimir forces 
for a system consisting of two spherical colloids immersed in a binary liquid mixture near its conso- 
lute point and facing a planar, homogeneous substrate. For several geometrical arrangements and 
boundary conditions we analyze the normal and the lateral critical Casimir forces acting on one of 
the two colloids. We find interesting features such as a change of sign of these forces upon varying 
either the position of one of the colloids or the temperature. By subtracting the pairwise forces from 
the total force we are able to determine the many-body forces acting on one of the colloids. We have 
found that the many-body contribution to the total critical Casimir force is more pronounced for 
small colloid-colloid and colloid-substrate distances, as well as for temperatures close to criticality, 
where the many-body contribution to the total force can reach up to 25%. 



I. INTRODUCTION 

The first step in describing the interaction between 
many particles is to determine their pair potential or 
the forces among a single pair. If the governing phys- 
ical equations are linear (like for gravity or electrostat- 
ics), this approach yields a quantitatively reliable de- 
scription of the physical system considered, based on 
the linear superposition principle. However, if nonlin- 
earities are present, linear superposition of pair poten- 
tials is no longer accurate and nonadditivity gives rise to 
many-body effects. These latter effects can lead, e.g., to 
a strengthening or weakening of the total force acting on 
a particle surrounded by more than a single other one, a 
change of sign of that force, or the appearance of stable 
or unstable configurations. Many-body effects appear in 
rather diverse systems such as nuclear matter [T], super- 
conductivity P], colloidal suspensions [31 2], quantum- 
electrodynamic Casimir forces [SUlOj. polymers [TT1 112j . 
nematic colloids [13] , and noble gases with van der Waals 
forces acting among them [T1HT7] . Each of these systems 
is characterized by a wide range of time and length scales. 
Integrating out the degrees of freedom associated with 
small scales (such as the solvent of colloidal solutes or 
polymers) for fixed configurations of the large particles, 
generates effective interactions among the latter, which 
are inherently not pairwise additive. This is the price 
to be paid for achieving a reduced description of a multi- 
component system. Driven by these effective interactions 
the large particles of the system may exhibit collective be- 
havior of their own (like aggregation or phase separation, 
see Refs. [TB] and [TH] and references therein), which can 
be described much easier if it is governed by pair poten- 
tials. In order to be able to judge whether this ansatz 
is adequate one has to check the relative magnitude of 
genuine many-body forces. 

In this paper we assess the quantitative influence 



of such many-body effects on critical Casimir forces 
(CCFs) [20-22 ]. These long-ranged forces arise as a con- 
sequence of the confinement of the order parameter fluc- 
tuations in a critical fluid |23j . They have been analyzed 
paradigmatically by studying the effective interaction be- 
tween a single colloidal particle and a homogeneous [2"4"1 - 
[3"2"| or inhomogeneous '■>'■'> 37 container wall as well as 
between two isolated colloidal particles [24-28] upon ap- 
proaching the critical point of the solvent. Here we add 
one sphere to the sphere-wall configuration, which is the 
simplest possibility to study many-body forces. (The wall 
mimics a third, very large sphere.) 

In order to be able to identify the latter ones one 
has to resort to a theoretical scheme which allows one 
to compute the forces between individual pairs and the 
three-body forces on the same footing. Since these forces 
are characterized by universal scaling functions, which 
depend on the various geometrical features of the con- 
figuration and on the thermodynamic state, we tackle 
this task by resorting to field-theoretic mean field theory 
(MFT), which captures the universal scaling functions as 
the leading contribution to their systematic expansion in 
terms of e = 4 — d spatial dimensions. Experience with 
corresponding previous studies for simple geometries tells 
that this approach does yield the relevant qualitative fea- 
tures of the actual universal scaling functions in d = 3; if 
suitably enhanced by renormalization group arguments 
these results reach a semi-quantitative status. We point 
out that even within this approximation the numerical 
implementation of this corresponding scheme poses a se- 
vere technical challenge. Thus at present this approach 
appears to be the only feasible one to explore the role of 
many-body critical Casimir forces within the full range 
of their scaling variables. 

Accordingly we consider the standard Landau- 
Ginzburg- Wilson Hamiltonian for critical phenomena of 
the Ising bulk universality class, which is given by 
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with suitable boundary conditions (BCs). In the case 
of a binary liquid mixture near its consolute (demixing) 
point, the order parameter cj>(r) is proportional to the 
deviation of the local concentration of one of the two 
species from the critical concentration. V is the volume 
accessible to the fluid, r is proportional to the reduced 
temperature t = (T — T c )/T c , and the coupling con- 
stant u > stabilizes the statistical weight exp(— Jff) 
in the two-phase region, i.e., for t < 0. Close to the 
bulk critical point T c the bulk correlation length £± di- 
verges as £±(f — > ± ) — £q \t\ ~ v , where v ~ 0.63 in 
d = 3 and v = 1/2 in d = 4, i.e., within MFT j3S]. The 
two non-universal amplitudes £ are of molecular size; 
they form the universal ratio £(j~/£o~ ~ 1-9 for d = 3 and 
£oV£o = \/2 for d = 4 [39]. The BCs reflect the generic 
adsorption preference of the confining surfaces for one of 
the two components of the mixture. For the critical ad- 
sorption fixed point [40] , the BC at each of the confining 
surfaces is either tf> — +oo or <fi = — oo, to which we refer 
as (+) or (— ), respectively. 

Within MFT the equilibrium order parameter distri- 
bution minimizes the Hamiltonian in Eq. ([I]) for the 
aforementioned BCs, i.e., Sjff [(/>]/ 5 (f> = 0. Far from any 
boundary the order parameter approaches its constant 
bulk value (<j>) = ±A \tf for t < or (0) = for t > 0. 
A is a non- universal bulk amplitude and j3 = 1/2 (for 
d = 4) is a standard critical exponent. In the following 
we consider the reduced order parameter m — u 1 ^ 2 (tp). 

The remainder of this paper is organized as follows. In 
Sec. [n] we define the system under consideration and the 
scaling functions for the CCFs as well as the normaliza- 
tion scheme. In Sec. |III| we present the numerical results 



where IZi < Ti-2 < • • ■ < are the semiaxes (or radii) 
of the hypercylinder and x = r\, y = r%, z = r^. If 
IZi = TZ% = ... = TZx> and d = D, the hypercylin- 
der reduces to a hypersphere. The generalization of d 
to values larger than 3 is introduced for technical rea- 
sons because d c = 4 is the upper critical dimension for 
the relevance of the fluctuations of the order parameter. 
These fluctuations lead to a behavior different from that 
obtained from the present MFT which (apart from loga- 
rithmic corrections [101 SI]) is valid in d = 4. We consider 
two hypercy finders in d = 4 with IZi — IZ2 = IZ3 = R 
and T> = 3. The two colloids are taken to be parallel 
along the fourth dimension with infinitely long hyper- 
axes in this direction. Considering hypercy finders, which 
are translationally invariant along the r^— axis, allows us 
to minimize j^[4>] numerically using a three dimensional 



obtained for the universal scaling functions of the CCFs, 
from which we extract and analyze the many-body ef- 
fects. In Sec. IIVI we summarize our results and draw 
some conclusions. 



II. PHYSICAL SYSTEM 

We study the normal and the lateral CCFs acting on 
two colloidal particles immersed in a near-critical binary 
liquid mixture and close to a homogeneous, planar sub- 
strate. We focus on the critical concentration which im- 
plies the absence of a bulk field conjugate to the order 
parameter [see Eq. Q]. The surfaces of the colloids and 
of the substrate are considered to exhibit a strong ad- 
sorption preference for one of the two components of the 
confined liquid leading to ± BCs. The forces are calcu- 
lated using the full three-dimensional numerical analysis 
of the appropriate MFT as given by Eq. Q. Specifi- 
cally, we consider two three-dimensional spheres of radii 
i?i and i?2 with BCs (a\) and (02), respectively, facing a 
homogeneous substrate with BC (b) at sphere-surface-to- 
substrate distances D\ and D2, respectively (see Fig. [I]). 
The coordinate system (cc, y, z) is chosen such that the 
centers of the spheres are located at (xi, 0, D\ +R\) and 
{x\ + Ri + L + R 2 , 0, D 2 + R2) so that the distance be- 
tween the centers, projected onto the cc-axis, is given by 
Ri+L+R 2 . The BCs of the whole system are represented 
by the set (01,02,6), where a±, 02, and b can be either 
+ or — . It is important to point out that we discuss 
colloidal particles with the shape of a hypercylinder 
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finite element method in order to obtain the spatially 
inhomogeneous order parameter profile 4>{x, y, z) for the 
geometries under consideration (see Fig. [T]) . 

In the case of an upper critical demixing point of the 
binary liquid mixture at the critical concentration, t > 
corresponds to the disordered (i.e., mixed) phase of the 
fluid, whereas t < corresponds to the ordered (i.e., 
phase separated) phase. The meaning of the sign is re- 
versed for a lower critical point. In the following we as- 
sume an upper critical point. 

The normal CCF i^fi 2 b) (D 1 ,D 2 , Ri,R 2 ,L, t) acting 
on sphere i in the presence of sphere j G {1-2} 

and i 7^ j) along the z-direction takes the scaling form 



C d ({U s }) = j r = (n,r a , . . . , r d ) G R d I J2 (^7) 
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FIG. 1: Two spherical colloidal particles of radii K\ and R% 
immersed in a near-critical binary liquid mixture (not shown) 
and close to a homogeneous, planar substrate at z = 0. The 
two colloidal particles with BCs (a\) and (02) are located at 
sphere-surface-to-substrate distances D\ and D2, respectively. 
The substrate exhibits BC (b). The lateral distance between 
the centers of the spheres along the :r-direction is given by 
R 1 + L + R2 , while the centers of both spheres lie in the plane 
y — 0. In the case of four spatial dimensions the figure shows 
a cut of the system, which is invariant along the fourth di- 
rection, i.e., the spheres correspond to parallel hypercylinders 
with one translationally invariant direction, which is Ta- 



where 5 = sign(i)L/£± (i.e., 3 = L/£ + for t > and 
S = — L/£- for t < 0). Note that the choice of 2 as the 
scaling variable does not depend on the type of particle 
the force acts on. Equation Q also describes the singular 
contribution to the lateral force near T c . The total CCF 
acting on particle i is 
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(oi,a 2 ,o) z 
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where and e z are the unit vectors pointing in x- and 
^-direction, respectively. Due to symmetry all other com- 
ponents of the CCF are zero. 

As a reference configuration we consider a single spher- 
ical colloid of radius R with BC (a) at a surface-to-surface 
distance D from a planar substrate with BC (b). This 
colloid experiences (only) a normal CCF 
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d-V+2 (ai,a 2 ,6) 



(e i ,A 1 ,A 2 ,n 1 ,n 2 ), (3) 



In the following we normalize the scaling functions 
and <;:i 2 ,6) by the amplitude K^Q = 
0, A = 1) of the CCF acting at T c on a single colloid 
for (+, +) BCs at a surface-to-surface distance D = R. 
Accordingly, in the following we consider the normalized 
scaling functions 



where Ax = D 1 /R 1 , A 2 = D 2 /R 2 , U ± = L/R u n 2 = 
L/R 2 , and 6, = sign(t)A/£± (i.e., 6, = A/£+ for <> 
and Qi = — A/£- for t < 0). Equation (|3| describes the 
singular contribution to the normal force emerging upon 
approaching T c . F^ 1 '^ is the force per length of the hy- 
percylinder due to its extension in the translationally in- 
variant direction. In the spirit of a systematic expansion 
in terms of e = 4 — d around the upper critical dimension 
we study the scaling functions K within MFT as given 
by Eq. ([T]) for hypercylinders in d = 4, which captures 
the correct scaling functions in d — 4 up to logarithmic 
corrections occurring in d = d c [40 1 41 , which we do not 
take into account here. Since MFT renders the leading 
contribution to an expansion around d = 4, geometrical 
configurations with small D/R, L/R, or -D/£+ are n °t ex ~ 
pected to be described reliably by the present approach 
due to the dimensional crossover in narrow slit-like re- 
gions, which is not captured by the e = 4 — d expansion. 

The colloidal particles will also experience a lateral 
CCF F { ( lf ao b) (D 1 ,D 2 ,R 1 ,R 2 ,L,t), for which it is con- 
venient to use the scaling form 
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with A s=z = ®i and A s=x = S. Experimentally it can be 
rather difficult to obtain Kfc*h(Q = 0, A = 1). A stan- 
dard alternative way to normalize is to take the more 
easily accessible amplitude A( + + ) for the CCF at T c be- 
tween two parallel plates with (+, +) BCs, which is given 
within MFT by (see Ref. [35] and references therein) 



A (+>+) - -24 



[/c(i/V§)] 



-283.6UT 
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where JC is the elliptic integral of the first kind. Within 
MFT the amplitude K ( ( * ,z ^(0 = 0, A = 1) can be ex- 
pressed in terms of A( + +y. 
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Equation ([£]) allows for a practical implementation of the 
aforementioned normalization, which eliminates the cou- 
pling constant u, which is unknown within MFT. 
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We calculate the normal and lateral forces directly 
from the numerically determined order parameter pro- 
files <j)(x, y, z) by using the stress tensor which, within 
the Ginzburg-Landau approach, is given by [210 H2] 



Due to the change of sign of K 
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with {k, 1} £ {x, y, z}. 

The first index of the stress tensor denotes the direc- 
tion of a force, the second index denotes the direction 
of the normal vector of the surface upon which the force 
acts. Therefore one has 
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where A; is a hypersurface enclosing particle i, hi is the 
i-th component (to be summed over) of its unit outward 
normal, and £d-v = J d d ~®r is the length of the d — 53- 
dimensional hyperaxis of C4. In particular we focus on 
the normal and lateral CCFs acting on colloid (2) for the 
configuration shown in Fig. [TJ for t ^ 0, and with the 
binary liquid mixture at its critical concentration. In the 
following analysis we consider colloid (1) to be fixed in 
space at a sphere-surface-to-substrate distance D\ = R\, 
equally sized colloids (i.e., R\ = R 2 ), and fixed (b = +) 
BC for the substrate. We proceed by varying either the 
vertical (z-direction) or the horizontal (x-direction) posi- 
tion of colloid (2) by varying either D 2 or L, respectively. 
We also consider different sets of BCs for the colloids. In 
the following results the numerical error is typically less 
than 5%, unless explicitly stated otherwise. 



III. RESULTS 

A. Total force 

In Fig. [2] we show the behavior of the normalized 
[Eq. 0] scaling function if(_' i£ ; 2i+) (e2,Ai = l,A 2 ,IIi = 
l,n 2 = 1) of the normal CCF acting on colloid (2) with 
(02 = ±) BCs close to a homogeneous substrate with 
(b = +) BC and in the presence of colloid (1) with 
(a 1 = — ) BC. The scaling functions are shown as func- 
tions of the scaling variable ratio O2/A2 = R 2 /£+, i.e., 
for t > 0. The various lines correspond to distinct fixed 
values of A 2 as the sphere-surface-to-substrate distance 
in units of the sphere radius. Thus Fig. [2] shows the tem- 
perature dependence of the normal CCF on colloid 2 for 
three different values of D 2 and for colloid (1) fixed in 
space. 

From Fig. [2] (a) one can see that, for colloid (1) with 
(01 = — ) BC and for any given value of <d 2 /A 2 , the 
scaling function of the normal CCF acting on colloid 
(2) with (a 2 = +) BC changes sign upon varying A 2 . 
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for any value of 



2 / A2 there is a certain value A 2 u; > 1 at which the 
normal CCF acting on colloid (2) vanishes. For sphere- 
surface-to-substrate distances sufficiently large such that 
A 2 > A^ ', colloid (2) is pushed away from the substrate 
due to the dominating repulsion between the two colloids 
in spite of the attraction by the substrate, whereas for 
A 2 < A 2 it is pulled to the substrate due to the dom- 
inating attraction between it and the substrate. This 
implies that, in the absence of additional forces, levita- 
tion of colloid (2) (i.e., zero total normal force) at height 
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A^R 2 is not stable against perturbations of the 
sphere-surface-to-substrate distance. On the other hand, 
upon varying temperature, any distance D 2 can become 
a stable levitation position for colloid (2) with (a 2 = — ) 
BC in the presence of colloid (1) with (ai = — ) BC [see 
Fig. [2] (b) , according to which each scaling function corre- 
sponding to a certain value of D 2 exhibits a zero so that, 
at this zero, increasing (decreasing) D 2 at fixed tempera- 
ture leads to an attraction (repulsion) to (from) the sub- 
strate]. In this case the attraction between the two col- 
loids is dominating for large sphere-surface-to-substrate 
distances A2, while the repulsion between colloid (2) and 
the substrate dominates for small values of A 2 . 

Figure [3] shows the behavior of the normalized scaling 
^ >+) (H,A 1 = l,A 2) n x = l,n 2 = l) of the 



(2,x) 
(ai,a 2 



function K 

lateral CCF acting on colloid (2) in the presence of col- 
loid (1) having the same BC, i.e., (ai) = (02) = (±)- The 
scaling functions are shown as functions of the scaling 
variable ratio S/n 2 = R 2 /£,+ - From Figs. [3] (a) and (b) 

(2 x ) 

one can infer that AT(±'± + j < 0. Therefore colloid (2) is 
always attracted towards colloid (1) which has the same 
BC. Hence the substrate does not change the sign of the 
lateral CCF as compared with the attractive lateral CCF 
in the absence of the confining substrate. However, the 
shapes of the scaling functions for (+, +, +) BCs [Fig. [3] 
(a)] differ from the ones for (-,-,+) BCs [Fig. |] (b)]; 
without the substrate, they are identical. In the former 
case and in contrast to the latter one, the scaling func- 
tions exhibit minima above T c , which is reminiscent of 
the shape of the corresponding scaling functions in the 
absence of the substrate. 

In Fig. [4] we show the results obtained for the normal- 

(2 x) 

ized scaling functions Krl + Ai = 1, A 2 = 1, Hi = 
n 2 ,n 2 ) of the lateral CCF acting on colloid (2). In 
Fig. (4^a) the scaling function is shown as function of the 
scaling variable ratio S/n 2 = R 2 /£,+ ; the black, red, and 
green curves correspond to L = R 2 , L = 1.5R 2 , and 
L = 2R 2 , respectively. In the absence of the substrate, 
the CCF between two colloids with opposite BCs is repul- 
sive. However, as shown in Fig. Qa), in the presence of 
the substrate with (b = +) BC, there is a change of sign in 
the scaling function of the lateral CCF. This implies that 
the lateral CCF acting between the two colloids changes 
from being attractive to being repulsive (or reverse) upon 
decreasing (increasing) the reduced temperature. Thus 
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2 /A 2 = R 2 /^ + 
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FIG. 2: Normalized scaling functions •K P (-',a 2 ,+) (02, Ai = 
1, A2,IIi = 1, II2 = 1) of the normal CCF acting on colloid 
(2) with BCs (02 = +) in (a) and (02 = -) in (b). The 
scaling functions are shown for t > as functions of the scaling 
variable ratio O2/A2 = i?2/£+ for three fixed values of the 
scaling variable A2 = -D2/-R2: A2 = 1 (black lines), 1.5 (red 
lines), and 2 (green lines), while Ai = IIi = II2 = 1 for all 
curves in (a) and (b) so that Di = Ri = R2 = L. For R2 
fixed the three curves correspond to three different vertical 
positions of colloid (2) with colloid (1) fixed in space. As 
expected, the forces become overall weaker upon increasing 

A2. -^(-^2,+) < (> 0) implies that the colloid is attracted 
to (repelled from) the substrate along the z-direction. 



temperature allows one to control both the strength and 
the sign of the lateral CCF in the case of two colloids 
with opposite BCs being near a wall. 

The at first sight unexpected lateral attraction between 
two colloids with opposite BCs in the presence of the sub- 

strate (i.e., K (_',+,+) < 0) can be understood as follows. 
In the absence of the two colloids, the order parameter 
profile m(r) is constant along any path within a plane 
z = const because in this case m(r) = m{z). Since the 
substrate area is much larger than the surface areas of the 
colloids, one can regard the immersion of these colloidal 
spheres as a perturbation of this profile. In Figs.^a) and 
(b) , the region within which the scaling function is nega- 
tive (corresponding to an attractive force) indicates that 
under these circumstances [i.e., when the colloids are suf- 
ficiently away from each other; see Fig. |4^b)] the pertur- 
bation generated by the presence of the spheres decreases 
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s/n 2 = R 2 /^ + 
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FIG. 3: Normalized scaling functions .KV,,' „ + \ (S, Ai = 
1, A2,IIi = l,n.2 = 1) of the lateral CCF acting on colloid 
(2) facing a homogeneous substrate with (6 = +) BC and in 
the presence of colloid (1) with (01) = (02) BCs. The scal- 
ing functions are shown for t > as functions of the scaling 
variable ratio E/II2 = i?2/£+ for three fixed values of the 
scaling variable A2 = D2/R2: A2 = 1 (black curves), 1.5 (red 
curves), and 2 (green curves), while Ai = IIi = II2 = 1 for 
all curves in (a) and (b) so that D\ = R% = R2 = L. For 
R2 fixed the three curves correspond to three different verti- 
cal positions of colloid (2) with colloid (1) fixed in space. As 
expected, the forces become overall weaker upon increasing 

A2. K (oi 01,+) < implies that colloid (2) is attracted to- 
wards colloid (l). Two different sets of (ai,a2,+) BCs are 
considered: (+,+,+) in (a) and (—,—,+) in (b). 



upon decreasing the lateral distance between them. This 
causes the colloids to move towards each other in order 
to weaken the perturbation by reducing its spatial exten- 

sion; this amounts to an attraction, i.e., K, |_ + j < 0. 

On the other hand, when they are sufficiently close to 
each other the pairwise interaction between the two col- 
loids dominates and the total lateral CCF is positive (i.e., 
repulsive). In Fig.[4|c) we show how the equilibrium lat- 
eral distance Lq measured in units of £ varies as function 
of temperature, i.e., £ for fixed R. 

In order to determine the preferred arrangement of the 
colloids, we have also analyzed the direction of the total 
CCF F?'^ )6) acting on colloid (2) [see Eq. ^} for sev- 
eral spatial configurations and BCs. For (+, +, +) BCs 
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FIG. 4: (a) Normalized scaling functions i£")J + Ai = 

1, A2 = 1, lli = 112,112) of the lateral CCF acting on colloid 
(2). Both colloids are taken to have the same size (ITi = II2 = 
II) and the same sphere-surface-to-substrate distance (Ai = 
A2 = 1). The scaling function is shown for t > as function 

(2 x) 

of the scaling variable ratio H/II2 = 7?2/£+- Kr_l^ t+ \ < 
(> 0) implies that colloid (2) is attracted (repelled) by col- 
loid (1). Black, red, and green curves correspond to lateral 
distances II2 = L/R2 = 1, 1.5, and 2, respectively, between 
the surfaces of the colloids [see Fig. [lj. The zero Hq(II2) of 



K 



(2,x) 
(-,+,+) 



implies that for any lateral distance L there is a re- 



duced temperature such that for £ + = Ho (II = LjR^L this 
distance represents an equilibrium lateral distance between 
the two colloids, provided they are located at equal sphere- 

— (2 x) 

surface-to-substrate distances, (b) K)_l as function of 
the reduced surface-to-surface distance II2 = L/R2 between 
the two colloids for E/II2 = R 2 /^+ = 3.162 (blue line), 3.873 
(yellow line), 4.472 (magenta line), and 5 (purple line). Each 
curve in (b) corresponds to the vertical dashed line with same 
color in panel (a). The change of sign of the scaling function 
(from positive to negative) as the distance between the col- 
loids is increased indicates that the lateral CCF changes from 
repulsive to attractive, which means that there is a lateral 
position Lq corresponding to a stable equilibrium point. In 
(c) we show how this equilibrium position, measured in units 
of £+, varies as function of temperature (i.e., £+) for fixed R 2 . 
Note that Lq is not proportional to £+, since £o/£+ is not 
a constant. As a guide to the eye the four data points are 
connected by straight lines. 




(1) 





(2) 



I + ± 



FIG. 5: Sketch showing the unit vector pointing into the di- 
rection of the resulting CCF acting on colloid (2) for (+, +, +) 
and (-,-,+) BCs, with Ai = IT = II2 = 1. The black 
rectangles represent the substrate and blue circles represent 
colloid (1) while black, red, and green circles represent colloid 
(2) with A2 = 1, 1.5, and 2, respectively. The centers of all 
colloids lie in the plane y = 0. 



we have found that the colloids tend to aggregate later- 
ally in such a way that several particles with (+) BC, 
facing a substrate with the same BC, can be expected 
to form a monolayer on the substrate. On the other 
hand, for the case of (— , — , +) BCs, we have found that 
the colloids can be expected to aggregate on top of each 
other so that a collection of colloids with such BCs is ex- 
pected to form three-dimensional sessile clusters. These 
tendencies become more pronounced upon approaching 
the critical point (see Fig. [5]) . Similar results have been 
found by Soyka et al. [33] in experiments using chemically 
patterned substrates. For them the authors have found 
indeed that colloids with (— ) BC distributed over those 
parts of the substrate with the same BC [which is equiva- 
lent to (+, +, +) BCs] aggregate and form a single layer. 
Moreover, they have found that colloids distributed over 
parts of the substrate with opposite BC [corresponding 
to (— . — , +) BCsl form three-dimensional clusters. 



B. Many-body forces 



We have determined the many-body force acting on 
particle (i) by subtracting from the total force F/*'*^ b ^ 
sec Eq. j5}] the sum of the pairwise forces acting 
on it, i.e., the colloid-colloid (CC) and the colloid- 
substrate (CS) forces. Accordingly the many-body CCF 
F {ai*a™b) {Di,D 2 ,Ri,R2,L,t) acting on colloid (i) is 
given by (see Fig. [T]) 



„(i,xz,MB) _ „(i,xz) , n n p R r s 



? (i,xz,CS) 
(Oi.6) 



(Di,Ri,t) . 



(12) 
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where 



F (aua2S^ Dl ' D2 ' Rl ' R2 ' L '^ = Tsign(aia 2 ) (e x cos a + e z sin a) x 

x f(a u a 2 ) (yJ{D 2 - Di + R~ 2 - Ri) 2 + (L + Ri + R 2 f -R x - R 2 ,R u R 2 ,t) (13) 



with 



a — arctan 



D 2 -D 1 +R 2 - Ri 
L + R 1 +R 2 



(14) 



is the pairwise colloid-colloid force (acting on colloid 2 
(— ) or 1 (+) with 2 having the larger ^-coordinate) ex- 
pressed in terms of the absolute value f( ai<a2 ) (u, i? 2 , t) 
of the force between two colloids at surface-to-surface dis- 



tance u in free space. F 



(i,xz,CS) 
(o»,6) 



(Di,Ri,t) is the CCF 



between the substrate and a single colloid (i). 
We have studied both the normal 



F, 



(i,z,MB) 



(01,02,6) 



and the lateral 



F, 



(i,x,MB) 
(ai,a 2 ,b) 



(ai,a 2 ,6) 



F (i,xz,MB) _ 
(ai,a 2 ,6) 1 



(15) 



(16) 



many-body CCFs which are characterized by correspond- 
ing scaling functions [compare Eqs. Q and (Eh]: 



p(i,z,MB) _ 
(ai,a 2 ,b) 



1?, 



(i,2,MB) ro 



d-S+2 (ai,a 2 ,fc) 



(e,,A 1 ,A 2 ,n 1 ,n 2 ), (17) 



and 



(i,x,MB) 
(ai,a 2 ,6) 



I?, 



^^'(S.A.A.n,,^). (18) 



Ld-Tl+2 (01,03,6) 



In Fig. [6] we show the normalized [see Eqs. ([6]) and Q] 
scaling functions K) ai +t+ \ (fc>2,Ai = 1,A 2 = l,lli = 
112,112) of the many-body normal CCF acting on colloid 
(2). This figure reveals similar results for (+, +, +) [Fig. [6] 
(a)] and (-,+,+) [Fig. |6] (b)] BCs. In these cases, each 
MB scaling function exhibits both a maximum and at 
least one minimum, the former one appearing for smaller 
values of the scaling variable R/£+ (i.e., at temperatures 
closer to T c ). For a certain range of temperatures close 
to T c , as the distance L between the colloids increases, 



the many-body normal CCF changes from attractive to 
repulsive. This shows that for each temperature within 
this range there is a lateral distance Lq IB /R for which 
the many-body contribution to the normal force acting 
on colloid (2) is zero. This means that under such condi- 
tions the sum of pairwise forces provides a quantitatively 
reliable description of the total force acting on colloid 
(2). For temperatures sufficiently far from T c , the many- 
body normal CCF is always attractive with a monotonic 
dependence on R/£+. Here, as in Figs.|]-EJ the CCFs 
decay exponentially for £ + — > 0. 

As expected, the many-body effects are more pro- 
nounced if the colloids are closer to each other and/or 
closer to the substrate. Indeed for situations in which 
the colloids are close to each other [see, e.g., the black 
curves m Figs. [6] (a) and (b)] we have found that when the 
normal many-body CCF reaches its maximal strength, 
corresponding to the minimum of the scaling function 

( 2 z Ad B ) 

' + + ^ at ~ 5, the relative contribution of the 

many-body CCF reaches 25% of the strength of the total 
normal CCF. For larger distances between the colloids 
[see, e.g., the green curves in Figs. [6] (a) and (b)] this rel- 
ative contribution is smaller (around 12% for II2 = 1.5 
and 6% for n 2 = 2). 

We are not aware of results for the quantum- 
electrodynamic Casimir interactions which are obtained 
along the same lines as our CCF analysis above. 
Nonetheless, in order to assess the significance of our 
results we compare them with the results in Ref. |10) . 
which is the closest comparable study which we have 
found in the literature. Therein the authors study theo- 
retically two dielectric spheres immersed in ethanol while 
facing a plate. Depending on the kind of fluid and on the 
materials of the spheres and of the plate as well as on 
the distances involved, also the quantum-electrodynamic 
Casimir force can be either attractive or repulsive. It 
is well known [331 |H] that for two parallel plates with 
permittivities e a and separated by a fluid with permit- 
tivity ef and without further boundaries, the quantum- 
electrodynamic Casimir force is repulsive if e a < < ei, 
within a suitable frequency range. In Ref. [5] it is stated 
that this also holds for two spheres immersed in a fluid. 
The authors of Ref. [10] analyze the effect of nonadditiv- 
ity for the above system by studying the influence of an 
additional, adjacent substrate on the equilibrium separa- 
tion d between two nanometer size dielectric spheres. To 



8 



normal MB force on (2) 




2 /A 2 = R 2 /^ + 

FIG. 6: Normalized scaling functions K^*^^\ (02, Ai = 
1, A2 = 1, III = n2,Il2) of the many-body normal CCF act- 
ing on colloid (2) for i?i = R 2 = R and Di = D 2 = R. The 
scaling functions are shown as functions of the scaling variable 
ratio O2/ A2 = R/i+ for the sets of BCs (+, +, +) in (a) and 
(—,+,+) in (b). The black, red, and green lines correspond 
to L/R — 1,1.5, and 2, respectively. Figures a) and |6|b) 
allow a direct comparison between the full CCF and the cor- 
responding many-body contribution (note the different scales 
of the ordinates.) 



this end, they consider two spheres of different materials 
with the same radii R\ = Rj = R and the same surface- 
to-plate distances Di = D 2 = D and analyze how the 
lateral equilibrium distance Lp between the spheres de- 
pends on D. By comparing the equilibrium distance Ljj 
with that in the absence of the substrate, L^, they find 
that Ld increases or decreases (depending on the kind of 
materials of the spheres) by as much as 15% as the dis- 
tance from the plate varies between D — 00 and D w R. 
They also find that "the sphere-plate interaction changes 
the sphere-sphere interaction with the same sign as D 
becomes smaller" , which means that if the sphere-plate 
force is repulsive (attractive), the sphere-sphere one will 
become more repulsive (attractive) upon decreasing the 
distance from the plate D. By construction these changes 
are genuine many-body contributions. In the case of two 
chemically different spheres, the sign of the many-body 
force contribution (i.e., whether it is attractive or repul- 
sive) agrees with the sign of the stronger one of the two 
individual sphere-plate interactions. 



In Fig. [7] we show schematically the system consid- 
ered in Ref. [10] [(a) and (b)] and the system considered 
here [(c), (d), and (e)]. For the quantum-electrodynamic 
Casimir effect, the dielectric spheres are represented by 
circles of equal radii, with the green one corresponding to 
a polystyrene sphere and the red one to a silicon sphere. 
The semi-infinite plates are represented by gray and yel- 
low rectangles for Teflon and gold, respectively. The 
whole configuration is immersed in ethanol which, for 
simplicity, is not shown in the figure. The dashed arrows 
indicate the direction of the strongest of the two pairwise 
sphere-substrate forces, while the solid arrows indicate 
the direction of the lateral many-body force. The direc- 
tions of the arrows in Figs. [7] (a) and (b) are chosen as 
to illustrate the findings in Ref. [TO], according to which 
the sign of the lateral many-body force is the same as the 
one of the strongest normal pairwise sphere-plate force: 
attractive in (a) and repulsive in (b). 

Also in the case of the critical Casimir forces, depicted 
in Figs. [7] (c), (d), and (e), we represent the colloids by 
circles and the laterally homogeneous semi-infinite sub- 
strate by rectangles. The orange filling represents the (+) 
BC while the blue filling represents the (— ) BC. Again, 
the dashed arrows indicate the direction of the stronger 
one of the two pairwise (normal) colloid-substrate forces, 
while the solid arrows indicate the direction of the lat- 
eral many-body contribution to the CCF. As one can 
infer from Fig. [8j the lateral many-body CCF acting on 
colloid (2) is always attractive for the given geometrical 
configuration, regardless of the BCs. 

We can also compare our results with those for two 
atoms close to the surface of a planar solid body. 
McLachlan [45] has tackled this problem by treating the 
solid as a uniform dielectric. By using the image method 
he derived an expression for the many-body corrections 
to the pairwise interaction energies, i.e., the atom-atom 
(London) and the atom-surface energies, in order to ob- 
tain the total interaction energy between the two atoms 
close to the surface. Qualitatively, he found that the 
leading contribution of the many-body correction leads 
to a repulsion if the atoms are side by side, i.e., at equal 
surface-to-substrate distances. 

Rauber et al. : 3B] used McLachlan's approach to study 
the electrodynamic screening of the van der Waals inter- 
action between adsorbed atoms and molecules and a sub- 
strate. The latter plays a role which is "analogous to that 
of the third body in the three-body interaction between 
two particles embedded in a three-dimensional medium" . 
The van der Waals interaction between the two atoms at 
equal distances from the substrate is altered by the pres- 
ence of the solid substrate and this perturbation is given 
by 051146] 

AT!"/ \ iCWl P\ C.S2 nQ , 

where p is the distance between the atoms, I is the height 
above the image plane, which is the same for both atoms, 
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quantum-electrodynamic Casimir interaction 



and p = 1 + 41 / p . The coefficients are given by 



(a) 



(b) 



0:0 



o 



critical Casimir interaction 

(c) (d) 







(e) 



o 



FIG. 7: Comparison of the direction of the lateral many- 
body forces (solid arrows) due to quantum-electrodynamic 
and critical Casimir interactions. In the case of the quantum- 
electrodynamic interaction [(a) and (b)] and alluding to the 
system studied in Ref. |10| . the circles represent the projec- 
tions of dielectric spheres with equal radii, the green one cor- 
responding to polystyrene and the red one to silicon; the rect- 
angles represent semi-infinite plates with their surfaces per- 
pendicular to the x — z plane, the gray and the yellow one 
being Teflon and gold, respectively. The system is immersed 
in ethanol, which is not indicated in the figure. In (a) and 
(b) each dashed arrow indicates the direction of the stronger 
one of the two corresponding pairwise forces between the di- 
electric spheres and the plate, which turns out to determine 
the direction of the many-body lateral force acting on the 
spheres: if the stronger one of the two pairwise forces is at- 
tractive [repulsive], the lateral many-body force will also be 
attractive [repulsive] (see Ref. [10]). Also in the case of the 
critical Casimir interaction [(c), (d), and (e)], the circles and 
rectangles represent projections of spherical colloids and of 
homogeneous substrates, respectively: orange and blue in- 
dicate (+) and (— ) BCs, respectively. In (c) and (d) the 
pairwise normal forces between each of the two spheres and 
the substrate are equal: attractive in (c) and repulsive in (d). 
In (e) the two pairwise normal forces have opposite directions 
with the repulsive one being the stronger one 32, 47, 48 . The 
corresponding dashed arrows have the same meaning as in (a) 
and (b). According to Fig.[8]the many-body lateral CCFs are 
attractive for all three cases(c), (d), and (e). The comparison 
shows that the systems in (a) and (c) behave similarly. How- 
ever, the behavior of system (b) has no counterpart for CCFs 
[see (d) and (e)]. In this figure all surface-to-surface distances 
equal the sphere radius, which in our notation corresponds to 
IT = n 2 = Ai = A 2 = 1. 



Csi 



:;/, 



a 2 (iOg(i()d( 



(20) 



and 



with 



ofc poo 

C S 2 = - / « 2 «)5 2 «R 

Jo 



<?«) = WO - 1] / WO + 1] , 



(21) 



(22) 



where lj — iQ is an imaginary frequency, a(oj) is the 
polarizability of the atoms (with the dimension of a vol- 
ume), and e(w) is the dielectric function of the solid (i.e., 
the substrate) . The lateral force due to the perturbation 
potential given by Eq. ( 19 1, which is the analogue of the 



many-body contribution to the lateral CCF, follows from 
differentiating AV with respect to p: 



dAV 
dp 



p7p3/2 



121 2 
VP 2 



2C. 



40l 4 C S i 
p 4 P 2 



si 



4C, 



CS2 
p372 



SI 



3Cg 2 

D 3/2 



(23) 



In Figure [9] we plot the lateral force given by Eq. ( 23 ) 
as function of the distance p between the two atoms for 
several (equal) distances I above the substrate. We use 
the values provided in Ref. [15] for the coefficients Csi 
and C S 2- Csi = 1.33 eVx(A) 6 and C 52 = 0.70 eVx(A) 6 , 
which correspond to Ne, and Csi = 17.70 eVx(A) 6 and 
Cs2 = 10.49 eVx(A) 6 , which correspond to Ar. As one 
can infer from Fig. [9j the many-body contribution to the 
lateral van der Waals force is always repulsive and, as the 
two atoms approach the substrate, its strength increases. 
On the other hand, in the case of the many-body contri- 
bution to the lateral CCF, we have found that it is at- 
tractive for all BCs considered, if the surface-to-surface 
distances between the spheres and the sphere-surface-to- 
substrate distances are equal to each other and to the 
radius of the spheres (see Fig.[8|. 

Further, we can quantitatively compare our results 
with those from Refs. [35J and [35] . To this end, we as- 
sign values to the geometrical parameters characterizing 
the configuration of the two atoms close to the substrate 
and compare the results of Refs. [45] and [46] with those 
for similar configurations in our model. For example, es- 
timating the many-body contribution to the lateral van 
der Waals force for a configuration of two atoms close to 
a substrate corresponding to the configuration associated 
with the black curve @5]. in Fig. [8] (i.e., L/R = D/R = 1 
in the case of the CCF and l/p = 1 in the case of the two 
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atoms), one obtains from Eq. (23) a value for the relative 
contribution of the many-body force to the lateral force 
which corresponds to ca. 15%. This is comparable with 
the relative contribution of the many-body force to the 
lateral CCF sufficiently close to T c , although in the case 
of the van der Waals force it is repulsive (Fig. [9]) while in 
the case of the CCF (Fig. [8| it is attractive. 

Considering the decay of the many-body contribution 
to the normal CCF as function of the surface-to-surface 
distance L between the spheres (for D% = D2 = D), 
we can compare the corresponding decay of the normal 
many-body force —dAV/dl given by the potential AV 
in Eq. (19). For small separations p between the atoms 



lateral MB force on (2) 



the many-body contribution to the normal van der Waals 
force increases as p -3 , while for large separations it de- 
cays as p~ s . By analyzing the data shown in Fig. [6] one 
finds that for 3 < i? 2 /C+ < 10 the scaling function of 
the many-body contribution to the normal CCF decays 
slower than L~ 2 . This means that in this temperature 
regime the many-body contribution to the normal CCF 
is much more long ranged than the corresponding contri- 
bution to the normal van der Waals force in the case of 
two atoms close to a surface. For fixed p the many-body 
contribution to the lateral van der Waals force decays 
as upon increasing the distance of both atoms from 
the substrate whereas the normal force on a single atom 
decays as l~ A . 

As a final remark we point out that we have not found 
a completely stable configuration for the two colloids 
(Fig. [T]) : whenever there is a stable position in the hor- 
izontal (vertical) direction, the force is nonzero in the 
vertical (horizontal) direction. For example, consider a 
vertical path with = 5 in Figs.[2jh) and[3jb). From 

the first one can see that along this path the normal 
CCF changes from being repulsive to being attractive as 
the sphere-surface-to-substrate distance for colloid (2) is 
increased, implying that there is a vertical position of 
colloid (2) in which the normal CCF is zero. However, 
according to Fig. [3](b) the lateral CCF is always attrac- 
tive regardless of the vertical position of colloid (2). This 
means that there is a configuration which is stable only 
in the normal direction. Accordingly, a dumbbell config- 
uration with a rigid thin fiber between the two colloids 
can levitate over the substrate. Whether this configura- 
tion is stable with respect to a vertical tilt remains as an 
open question. 



0- 




A l = A 2 = n i= n 2 =1 



_l I I I I I I I L 



10 

s/n 2 = R 2 /^ + 



FIG. 8: (a) Normalized scaling functions K^^'^^l(S, Ai = 
1, A2 = 1, IIi = 1, II2 = 1) of the lateral many-body CCF act- 
ing on colloid (2) for Ri = R2 — R and Di — D2 = L = R. 
The scaling function is shown as function of the scaling vari- 
able ratio H/II2 = R2H+ — R/£+- The black, red, and 
green lines correspond to the BCs ( — ,—,+), (—,+,+), and 
(+, +, +), respectively. For all three BCs the many-body con- 
tribution is not monotonic as function of temperature. Quan- 
titatively the green, red, and black curves here should be com- 
pared with the black curves in Figs. |3|a), |3|b), and[4ja), 
respectively. However, for the data shown in this figure the 
error bars (not shown) due to limits of the numerical accuracy 
are between 10% and 15%. This is the main reason why we 
refrain from showing what would be an instructive plot such 

as if(a^^+) as a function of II = IT = n2 for various values 
of A == Ai == A2 and £ (as we did in Fig. B, which would 
allow for a direct comparison with the case of atoms. 
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C. Vertical alignment 

We now turn our attention to the case in which the 
colloids are vertically aligned with respect to a planar, 
homogeneous substrate, i.e., when their centers have the 
same coordinates in both the x and the y direction (see 
Fig. 10). We focus on the normal CCF acting on colloid 



(1) when the system is immersed in a near-critical binary 
liquid mixture at its critical concentration. As before 
we consider ± BCs corresponding to a strong adsorption 



FIG. 9: Many-body contribution to the lateral van der Waals 
force between two atoms adsorbed at a dielectric, homoge- 
neous substrate [Eq. ( |23[ )], derived from the expression for 
the excess potential given in Refs. [45] and [46] [Eq. ( 19 1] . 
The forces are plotted as functions of the lateral separation p 
between the two atoms for several equal vertical distances I of 
the atoms from the substrate. The curves correspond to two 
sets of values for the coefficients Csi and Cs2 in Eq. ( |19[ ) 46 , 
corresponding to Ne (solid lines) and Ar (dashed lines). 
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FIG. 10: Two spherical colloidal particles of radii Ri and R2 
immersed in a near-critical binary liquid mixture (not shown) 
and close to a homogeneous, planar substrate at 2 = 0. The 
colloidal particle (1) with BC (ai) is located vertically at the 
sphere-surface-to-substrate distance D\ , whereas the colloidal 
particle (2) with BC (02) is located vertically at the surface- 
to-surface distance S between the spheres. The substrate ex- 
hibits BC (6). The vertical distance between the centers of the 
spheres along the z-direction is given by R1 + S+R2, while the 
centers of both spheres lie on the vertical axis x — = y. In 
the case of four spatial dimensions the figure shows a three- 
dimensional cut of the system, which is invariant along the 
fourth direction, i.e., the spheres correspond to parallel hyper- 
cylinders with one translationally invariant direction, which 
is r 4 [see Eq. p|] . 



preference for one of the two components of the confined 
liquid. In particular, we consider two three-dimensional 
spheres of radii R\ and R2 with BCs (ai) and (02), re- 
spectively, facing a laterally homogeneous substrate with 
BC (b). Colloid (1) is positioned at a sphere-surface-to- 
substrate distance D\ and colloid (2) is at a surface-to- 
surface distance S from colloid (1) (see Fig. [To] ). The 
coordinate system (x,y,z) is chosen such that the cen- 
ters of the spheres are located at (0,0, D\ + Ri) and 
(0, 0, Di + 2Rx + S + R2) so that the distance between 
the centers, along the z-axis, is given by R\ + S + R 2 . As 
before, the BCs of the system as a whole are represented 
by the set (a\, a 2 , b), where ai, a 2 , and b can be either + 
or — . 



The normal CCF F, 



(a 1 ,a 2 ,b) 



(Di, S, Ri, R 2 , t) acting on 



colloid (1) along the ^-direction takes the scaling form 



F,^' z _ L JDi,S,Ri,R2,t) 



(ai,a 2 ,6) 1 



Ri 



where a = sign(i)5/£± (i.e., er = S/£ + for t > and 
a = -S/Z- for t < 0), A x = Di/R u A 2 = D x /R 2 , 



and Si = S/Ri; S 2 = S/R 2 = S1A2/A1. Equation ((24j) 
describes the singular contribution to the normal force 
emerging upon approaching T c . F 1 -- 1 '^ is the force on a 



hypercylinder divided by its extension in the translation- 
ally invariant direction [see Eq. pj)]. We use the same 
reference system as the one described by Eq. ([6| in or- 



der to normalize the scaling function defined in Eq. (24) 
according to 



K 



(M) 

(ai,a 2 ,fc) 



(a,A 1 ,A 2 ,<5i) 



K, 



(♦>*) 

(+,+: 



(e = o,A = i) 



(25) 



We calculate 
to the normal 



F, 



(l,z,MB). 



the many-body contribution 

CCF acting on particle (1) 

(ai ai b) C^l> Rii -^2) t) by subtracting from the 

total force the sum of the pairwise forces acting on it, 
i.e., the colloid-colloid and the colloid-substrate forces 
[see Eq. 
form 



This many-body force takes the scaling 



F, 



(1,2, MB) 



(a 1 ,a 2 ,b) 



k R T 



S d-V+2 



r(l,z,MB)t 
l (ai,a 2 ,b) 



(26) 



In Fig. 11 we show the normalized [see Eqs. j§, 
(25), and (26)] scaling functions K^'^^(a, Ai, A 2 = 
Ai,<5i = 1) of the many-body contribution to the nor- 
mal CCF acting on colloid (1) as functions of the scaling 
variable ratio a/81 = Ri/£+ f° r t w0 spherical colloids of 
the same size (R\ = R 2 = R). Keeping the surface-to- 
surface distance between the spheres fixed at S — R, we 
vary the sphere-surface-to-substrate distance Di for sev- 
eral BCs: (+, +, +) in (a), (+, +) in (b), (-, +, +) in 
(c), and (-,-,+) in (d). 

From Figs. [Tl|a) and (d) one can infer that if the col- 
loids have symmetric BCs, the scaling function of the 
many-body normal CCF acting on colloid (1) is negative 
(i.e., it is directed towards the substrate) for any value 
of R/£+. On the other hand, for non-symmetric BCs be- 
tween the colloids [Figs. [lT|(b) and (c)], the many-body 
contribution to the normal CCF acting on colloid (1) is 
positive for any value of -/?/£+• The apparent change of 
sign in Figs.[TT|b) and (c) is likely to be an artifact occur- 
ring within the error bars due to numerical imprecision. 
The relative contribution of the many-body CCF to the 
total force is between 10% and 15%. 

We point out that this configuration with the the two 
colloids vertically aligned with respect to the substrate 
allows for a wide range of interesting aspects which will 
be further explored in future works. 



IV. CONCLUSIONS AND DISCUSSION 

We have investigated critical Casimir forces (CCFs) for 
a system composed of two equally sized spherical colloids 
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(i?i = i?2 = R) immersed in a near-critical binary liquid 
mixture and close to a laterally homogeneous substrate 
with (b = +) boundary condition (BC) (see Fig. [lj. By 
denoting the set of BCs of the system as (oj, 0,2, b), where 
ai corresponds to the BC at colloid i = 1, 2 and b to the 
BC at the substrate, we have first focused on the total 
normal and lateral forces acting on one of the colloids 
[labeled as "colloid (2)"] for several geometrical configu- 
rations of the system and various combinations of BCs 
at the colloids. Both the normal and the lateral forces 
are characterized by universal scaling functions [Eqs. ^ 
and Q, respectively], which have been studied in the 
one-phase region of the solvent as functions of R2/C+ 
and L/R2. L is the surface-to-surface distance between 
the two colloids, and is the bulk correlation length of 
the binary mixture in the mixed phase. We have used 
mean-field theory together with a finite element method 
in order to calculate the order parameter profiles, from 
which the stress tensor renders the normalized scaling 
functions associated with the CCFs. 

For the scaling function of the total normal CCF act- 
ing on colloid (2) with (0,2 = +) BC, in the presence 
of colloid (1) with (ai = — ), we have found (Fig. [2]ja)) 
that the scaling function changes sign for a fixed value 
of i?2/^+ as the distance D2 between colloid (2) and the 
substrate increases, signaling the occurrence of an unsta- 
ble mechanical equilibrium configuration of a vanishing 
normal force. For the total normal CCF acting on col- 



loid (2) with (0,2 = — ) BC and in the presence of colloid 
(1) with (ai = — ), we have found (Fig. ^b)) that the 
force changes sign upon changing the temperature. For 
this combination of BCs, the equilibrium configuration 
of colloid (2) is stable in the normal direction. 

Without a substrate, at the critical composition of the 
solvent the CCF between two (+) spheres is identical to 
the one between two (— ) ones. This degeneracy is lifted 
by the presence of the substrate as one can infer from the 
comparison of the scaling functions for the lateral CCFs 
for (+,+,+) and (-,-,+) BCs [see Figs. a) and (b), 
respectively]. In the first case, the shape of the scaling 
function resembles that of the two colloids far away from 
the substrate, with a minimum at i?2/£+ ~ 2.5. In the 
second case this minimum at T > T c disappears. These 
substrate-induced changes are more pronounced if the 
two spheres are close to the substrate (Fig. |3| . Without 
a substrate the CCF between spheres of opposite BCs 
is purely repulsive. In the presence of a substrate the 
corresponding lateral CCF for (—,+,+) BCs turns at- 
tractive for large lateral distances L (Fig. which is a 
pure many-body effect. 

We have also studied the direction of the total CCF 
acting on colloid (2) for various spatial configurations 
and BCs in order to assess the preferred arrangement of 
the colloids. For (+, +, +) BCs we have found that they 
tend to aggregate laterally. In this case a collection of col- 
loids with (+) BCs, facing a substrate with the same BC, 
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are expected to form a monolayer on the substrate. For 
the situation of (—,—,+) BCs, we have found that the 
colloids are expected to aggregate on top of each other. 
This indicates that a set of several colloids with such BCs 
are expected to form three-dimensional clusters. These 
tendencies are enhanced upon approaching the critical 
point (Fig. 5]). 

By calculating the pairwise colloid-colloid and colloid- 
substrate forces and subtracting them from the total 
force, we have determined the pure many-body contri- 
bution to the force acting on colloid (2). For the scaling 
functions associated with the normal many-body CCFs 
we have found the interesting feature of a change of sign 
at fixed temperature upon varying the lateral position 
of colloid (2) (Fig. |6|. This implies that, for a given 
temperature, there is a lateral position where the normal 
many-body CCF is zero, in which case the sum of pair- 
wise forces provides a quantitatively reliable description 
of the interactions of the system. As expected we have 
found that the contribution of the many-body CCFs to 
the total force is large if the colloid-colloid and colloid- 
substrate distances are small, as well as if the binary 
liquid mixture is close to its critical point. 

We have compared our results with corresponding ones 
for quantum-electrodynamic Casimir interactions. To 
this end we have referred to the results in Ref. [TU] for 
two dielectric spheres immersed in ethanol and facing a 
plate. These authors analyze the influence of the dis- 
tance D = D\ = D2 from the plate on the equilibrium 
separation Ld between the spheres, which are subject to 
quantum-electrodynamic Casimir forces. They find that 
the lateral many-body force is attractive (repulsive) if 
the stronger one of the two normal sphere-plate forces is 
attractive (repulsive) [Figs. [7^a) and (b)]. On the other 
hand, in the case of CCFs we have found that for a config- 
uration in which the surface-to-surface distance between 
the colloids is equal to the sphere-surface-to-substrate 
ones and equal to the radius of the spheres, the many- 
body contribution to the lateral CCF is always attractive, 
regardless of the BCs (Fig. [8]). 

We have also compared our results with the corre- 



sponding ones for the case of two atoms close to the 
planar surface of a solid body. In this respect we have 
referred to the McLachlan model [IS] for the many-body 
contribution to the van der Waals potential [see Eq. ( 19 ) 
and Fig. [9j and the results from Ref. jH] . From this com- 
parison we have found that if the two atoms are fixed at 
the same distance from the surface of the solid body, 
the normal many-body contribution to the total van der 
Waals force decays with the atom-atom distance p as p~ 8 
for large atom-atom distances. This decay is much faster 
than the decay we estimate for the many-body contribu- 
tion to the normal CCF, which within a suitable range 
appears to be slower than L~ 2 . Furthermore, we have 
found that the many-body contribution to the lateral 
van der Waals force is repulsive while the correspond- 
ing many-body CCF is attractive regardless of the set of 
BCs. 

Finally we have considered the configuration in which 
the two colloids are vertically aligned with respect to the 
substrate (Fig. 10). We have calculated the many-body 
contribution to the normal CCF acting on colloid (1) 
for two spherical colloids of the same size (Ri = R2 = 
R) keeping the sphere-surface-to-surface distance S = R 
fixed (Fig. 11). We have varied the sphere-surface-to- 



substrate distance D\ for several BCs and have found 
that if the colloids have the same BCs, the many-body 
contribution to the normal CCF is directed towards the 
substrate [Figs. [Tl|a) and (d)], whereas for colloids with 
opposite BCs, the many-body contribution to the normal 
CCF is directed away from the substrate [Figs.[TT|b) and 
(c)]. We have found that the contribution of the many- 
body CCF to the total force is between 10% and 15%. 
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